DFTT 57/95 
September 1995 



c3 



X 



Recent results in high temperature Q 
Lattice Gauge Theories 



M. Caselle 



, Dipartimento di Fisica Teorica dell'Universita di Torino 

Istituto Nazionale di Fisica Nucleare, Sezione di Torino 
via P.Giuria 1, 1-10125 Torino, Italy f\ 

On 

o i 

Abstract 

o 

\^ . We review some analytic results on the deconfinement transition in pure 

lattice gauge theories. In particular we discuss the relationship between the 
deconfinement transition in the {d + l)-dimensional SU(2) model and the 
magnetization transition in the d-dimensional Ising model. This analysis 
leads to a precise estimate of the deconfinement temperature which agrees 
well with that obtained with a Montecarlo simulation in the case in which 
the lattice has only one link in the compactified time direction. 
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1 Introduction 



One of the most interesting predictions of QCD is the existence of a deconfinement 
transition at some critical temperature T c . Finding a precise characterization of 
this transition is, however, still an open problem. In particular there are two main 
open questions: the first one is the identification of the order of the transition 
and, in case it is of second order, of its critical indices. The second one is the 
precise location of the deconfinement temperature. The natural framework to pose 
these questions is that of the finite temperature Lattice Gauge Theories (LGT). In 
this framework, a seminal contribution was given more than ten years ago by B. 
Svetitsky and L.Yaffe |]J in the case of pure gauge theories. They showed that, if 
the deconfining transition of a given (d + 1) dimensional gauge theory is of second 
order, then its universality class should coincide with that of the d- dimensional spin 
model , with symmetry group the center of the original gauge group. This result is 
usually known as the "Svetitsky- Yaffe (SY) conjecture" and has been confirmed in 
these last ten years by several Montecarlo simulations. It must be noticed however 
that the SY conjecture gives no information on the location of the deconfinement 
transition. Trying to answer to this last question will be the main goal of the present 
contribution. In particular we shall see how far one can go in trying to estimate 
the critical temperature by using only analytic methods. We shall concentrate only 
on pure gauge theories with gauge group SU(2), but most of our results can be 
straightforwardly extended to SU(N) models with N > 2. 

During these last years the best estimates of the critical temperature have been 
obtained by means of Montecarlo simulations, which are certainly the most powerful 
tool to extract quantitative results from LGT. However we think that it is important 
in itself to have some independent analytic estimate of the location of the critical 
point, besides the outputs of the computer simulations, to reach a deeper theoretical 
understanding of the deconfinement transition. The attempts to obtain analytically 
the critical temperature have a rather long history, starting more than ten years 
ago 0, D, H, D, §]. However the strategy has always been essentially the same: first 
construct an effective action in terms of the Polyakov loops (which, as we shall 
see below, are the relevant dynamical variables in the physics of deconfinement for 
pure gauge theories). Second, use a mean field approximation to extract the critical 
coupling. A common feature of all these attempts was that the effective actions were 
always constructed neglecting the spacelike part of the action. As a consequence it 
was impossible to reach a consistent continuum limit for the critical temperature. 
Moreover, as a consequence of the mean field approximation, the estimates of the 
critical temperature were in general affected by large systematic errors. 

The aim of the present contribution is to show that it is indeed possible to over- 
come these two difficulties. First we shall construct in the SU(2) case an improved 
effective action which takes also into account the spacelike part of the original Wil- 
son action. Second we shall avoid the mean field approximation, and shall instead 
obtain the critical temperature by mapping (following the SY conjecture) the gauge 
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theory to a suitable Ising-like model, and then using the fact that the critical tem- 
perature of the Ising model is known exactly in d = 2, and with very high precision 
in d = 3. Let us stress that this is not the only possible way to avoid the mean field 
approximation. Another interesting possibility is to study the SU(N) models in 
the N —>■ oo limit where one can use some results recently obtained in the context 
of random matrix models and two-dimensional exactly solvable gauge theories to 
discuss the deconfinement transition and fix the critical temperature. We shall not 
describe here this approach, the interested reader can find in a comprehensive 
discussion on the subject. 

This contribution is organized as follows. Sect. 2 will be devoted to a brief 
introduction to finite temperature LGT and to the Svetitsky- Yaffe conjecture. In 
sect. 3 we shall construct the effective action, which we shall then use in sect. 4 to 
extract the critical deconfinement temperature. Finally sect. 5 will be devoted to 
some concluding remarks. 



2 Finite Temperature LGT 
2.1 General setting 

Let us consider a pure gauge theory with gauge group SU(N), defined on a d + 1 
dimensional cubic lattice. In order to describe a finite temperature LGT, we have to 
impose periodic boundary conditions in one direction (which we shall call from now 
on "time-like" direction), while the boundary conditions in the other d direction 
(which we shall call "space-like") can be chosen freely. We take a lattice of N t (N s ) 
spacings in the time (space) direction, and we work with the pure gauge theory, 
containing only gauge fields described by the link variables U n $ G SU(N), where 
n = (x, t) denotes the space-time position of the link and i its direction. It is useful 
to choose different bare couplings in the time and space directions. Let us call them 
fit and fi s respectively. The Wilson action is then 

S W = T, 4 Re I A E Tr f (Un;0i) +&E Tr f (Un;ij) \ , (1) 

where Trf denotes the trace in the fundamental representation and U n -oi (U n -ij) are 
the time-like (space-like) plaquette variables, defined as usual by 

U n -ij = U n - t iU n +i-jU n +j.jU n .j . (2) 

In the following we shall call S s (St) the space-like (time-like) part of Sw- fis 
and fit are related to the (bare) gauge coupling g and to the temperature T by the 
usual relations 



where a is the space-like lattice spacing, while is the time-like spacing. 

In a finite temperature discretization it is possible to define gauge invariant 
observables which are topologically non-trivial, as a consequence of the periodic 
boundary conditions in the time directions. The simplest choice is the Polyakov 
loop, defined in terms of link variables as 

N t 

P(£)=Tr { P 3 = Tr { ]l(y s ,w) . (4) 
t=i 

In the following we shall call Pg, "open Polyakov line". 

As it is well known, the finite temperature theory has a new global symmetry 
(unrelated to the gauge symmetry), with symmetry group the center C of the gauge 
group (in our case Z 2 ). The Polyakov loop is a natural order parameter for this 
symmetry. 

In d > 1, finite temperature gauge theories admit a deconfinement transition 
at T = T c , separating the high temperature, deconfined, phase (T > T c ) from 
the low temperature, confining domain (T < T c ). The high temperature regime is 
characterized by the breaking of the global symmetry with respect to the center of 
the group. In this phase the Polyakov loop has a non-zero expectation value, and 
it is an element of the center of the gauge group. In the low temperature phase 
the center symmetry is conserved and the expectation value of the Polyakov loop 
is zero. The relevant feature of the Polyakov loop is that at the same time it is also 
the order parameter for the deconfinement transition. In fact its expectation value 
is related to the free energy of an isolated, static quark as follows: 

< P > oc exp(-F 3 ) . (5) 

As a consequence, in the low temperature phase it would require an infinite energy to 
create from the vacuum an isolated quark. Hence in this phase quarks are confined. 
On the contrary, in the high temperature phase isolated quark can exist: this is 
the deconfined phase. The critical point in which the center symmetry is broken 
can thus be interpreted as the deconfinement transition. The corresponding critical 
temperature T c will be denoted in the following as the deconfinement temperature. 

2.2 Svetitsky-Yaffe conjecture 

The idea on which the SY conjecture is based is that, if one would be able to in- 
tegrate out all the gauge degrees of freedom of the original (d + l)-dimensional 
model except those related to the Polyakov loops then the resulting effective the- 
ory for the Polyakov loops would be a d- dimensional spin system with symmetry 
group C. The deconfinement transition of the original model would become the 
order-disorder transition of the effective spin system. This effective theory would 
obviously have very complicated interactions, but Svetitsky and Yaffe were able to 
argue that all these interactions should be short ranged. As a consequence, if the 
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transition point of the effective spin system is of second order, near this critical 
point, where the correlation length becomes infinite, these short ranged interac- 
tions can be neglected, and the universality class of the deconfinement transition 
should coincide with that of the simple spin model with only nearest neighbour 
interactions and the same global symmetry group. As a consequence all the criti- 
cal indices describing the two transitions and all the adimensional ratios of scaling 
quantities should coincide in the limit. In particular the (d+ l)-dimensional SU(2) 
gauge theory, which is known to have a second order deconfinement transition, is 
characterized by the same renormalization group fixed point of the <i-dimensional 
Ising model. 

Unfortunately the SY argument alone cannot help to fix the critical temperature, 
since the precise mapping between the Ising coupling and that of the original gauge 
model, requires taking into account exactly those short ranged interaction which 
we neglected above. 



2.3 The SU (2) case: character expansion 

In the following we shall concentrate on the SU (2) model. There are two important 
features which greatly simplify the analysis in this case. The first one is that 
according to the SY conjecture the model can be mapped, at the deconfinement 
point into the spin Ising model, which is exactly solved in d = 2 and very well known 
in d — 3. The second important feature is that in the £77(2) case the character 
expansion (which plays an important role in the construction of the effective action) 
is very easy to handle. 

Let us briefly summarize few results. The character of the group element U in 
the j th representation is: 

where Trj denotes the trace in the j th representation and 9 is defined according to 
the following parametrization of U in the fundamental representation: 

U = cos(0)l + iansin(6) (7) 

where n is a tridimensional unit vector and <7j are the three Pauli matrices. Notice, 
as a side remark, that with this parametrization the Haar measure has the following 
form: 

DU = 8 m 2 (e)-^- (8) 

and the Polyakov loop becomes P(x) = 2 cos(6g) 

The following orthogonality relations between characters hold: 

J D U X r(U) X s(U) = 5r,s (9) 
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Y,drXr(U V" 1 ) = 8{U,V) (10) 

r 

where d r denotes the dimensions of the r th representation: d r — 2r + 1. In the 
following we shall use two important properties of the characters: 

J D U X r(U) XsiU^V) = Sr,s*^p- , (11) 

/ D UxAUVxU-^) = ^-Xr(Vx)xr(V 2 ) ■ (12) 
J d r 

The character expansion of the Wilson action has a particularly simple form: 
e f^) = ^ 2 (2i + l)^® X ,(C7), j = 0,1,1... (13) 

where I n (/3) is the n th modified Bessel function. 

In the following we shall often use the normalized version of the character ex- 
pansion in which the coefficient of the trivial representation is set to 1. 

el T ^) = G(/5)E(2j + l)^^X,(f/), J = 0,i,l... (14) 

where G{(3) = 2Ji(/3)//3 is an irrelevant constant that we shall often omit. 

3 Construction of the Effective Action 

In this section our goal is to construct an effective action for the finite temperature 
LGT in terms of the Polyakov loops only. This implies that we must be able to 
integrate out exactly all the spacelike variables so that the only remaining degrees 
of freedom at the end are exactly the Polyakov loops. Notice that in this way 
the resulting effective action would live in d dimensions (one dimension less than 
the starting model). This is exactly along the line of the original Svetitsky-Yaffe 
program. In trying to follow such a program one must necessarily make some 
approximation. In order to obtain a good approximation of the original Wilson 
action, one must identify the physically relevant degrees of freedom, and then try to 
keep them unchanged when constructing the effective action. Following [p]] and 0, 
we assume that the physics of the deconfinement transition is dominated by the 
timelike plaquettes, and try to keep as far as possible unchanged this part of the 
original action. Accordingly we treat the spacelike part of the Wilson action: S s as a 
perturbation of the timelike part St and take care of the contributions coming from 
S s by making a strong coupling expansion in f3 s . The main difference with respect 
to the usual approach is that in this case the time-like part of action is treated 
exactly or, equivalently, that the expansion in (3 t is summed up to all orders. The 



5 



only remaining expansion parameter is thus (3 S . In particular the zeroth order in (3 S 
will contain the timelike plaquettes only. It is not at all obvious that the integration 
over the spacelike links could be done to all orders in fit, but it turns out that it 
can be done exactly in the framework of the characters expansion order by order 
in (3 S . In particular we shall discuss the zeroth and first order in (3 S only, which will 
be enough to our purposes, but there is in principle no obstruction to go to higher 
orders. The result for any given order in (3 S can be expressed as an infinite sum 
over characters. 

Remarkably enough in the N t = 1 case this series can be summed exactly and 
the result can be written in a closed form. This is essentially due to the fact the 
if Nf = 1 this same effective action can be obtained in a completely different way, 
using techniques typical of matrix models (see below), thus allowing a non trivial 
check of all our strong coupling results. Another interesting feature of the N t = 1 
limit is that in this case very precise Montecarlo estimates of T c exist, with which we 
can compare our analytic predictions. For instance, in (3+1) dimensions the critical 
coupling (3 C at which deconfinement occurs is estimated to be: (for the SU(2) model 
with N t — 1) (3 C — 0.8730(2) ||. Such an impressive precision is due to that fact 
that in the N t = 1 case (and only in this case) one can simulate the gauge model by 
using a cluster non-local algorithm (see || for the details). This makes the SU(2), 
Nt = 1 model a perfect laboratory to test our techniques, and we shall concentrate 
on this case in sect. 4 . We shall briefly comment on the extension of our results to 
N t > 1 in sect. 5 . 



3.1 Expansion in (3 S of the effective action 



The effective action S e s for the Polyakov loops Pg = Titti Vx is obtained integrating 
over all the spacelike degrees of freedom in the action (fj). As explained previously, 
our approach is to consider the contributions from the spacelike plaquettes up to 
a certain order in f3 s only. So, for our purposes, it will be convenient to expand 
separately the spacelike and the timelike part of the action ([!]): 



J ]J DU x,t;i exp S w 

x,t;i 

n n 



x,t:i 



x,t:i 



1 + \- d.hi+l^Ay.m ^ 



>x±Mi 



X 



x,t:Li \ 7=1 M/VJ 



(15) 



Specifically, we work out here the effective action up to 0(f3l). This means that 
in eq. (|i~5|) we must look only at the terms containing at most a single space- 
like plaquette in the adjoint representation, Xi{Ux,t;ij) , or two space-like plaquettes 
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in the fundamental, x±(Ux,ti;ij)xi (Uff,t 2 ;ki)- Due to the orthogonality relations for 
characters, it easy to convince oneself that a pair of plaquettes in the fundamental 
representation do actually contribute to the integral only if they appear in the 
same spatial position (at two different times ti and t 2 ); for the same reason a 
single fundamental plaquette cannot contribute. We are thus lead to the following 
expression: 



exp(S cff ) = J J] DU m Jl 1 + E dj-f 

x,t;i x,t;i \ j = - 



2j+l {^t 



(A 



(fit 



t:0i) 



x 11+ E 

x,i<j 



E3^ 

t=l J l 



In the following we shall calculate these integrals for a generic value of N t 



if 

*Kt2 V J 1 



Xl(^?,t 1 rfi)Xi(^,teM)( 1 p) 



3.2 Zeroth order approximation 

Let us first study the contribution which gives the 0(/3°) result and corresponds to 
the "1" in the second factor of eq.([H]). In this case the integral only contains the 
timelike part of the Wilson action: 



exp(S'o) 



n 



, (17) 



and it is easy to integrate all the spacelike links. The reason is that each spacelike 
link only belongs to two timelike plaquettes, hence by making a character expansion, 
it can be exactly integrated out. Let us do this integration in two steps, for future 
commodity. First let us integrate (by using eq.(]TT|)) all the spacelike links except 
the lowermost ones (which, due to the periodic boundary conditions coincides with 
the uppermost). We obtain: 



exp(S )=n 1 + E 



(Uz-iPg+iUliPt 



\ 



J 



;is) 



where Pg is the open Polyakov line (whose trace is the Polyakov loop) in the site x 
and Ug-i are the remaining lowermost spacelike links. Integrating also on Ug-i (this 
time, by using eq.flTTD) we end up with 



exp(S'o 



^ oo 

n !+e 



hj+i (fit 
h(Pt) 



Xj (p s+i ) xA4) 



\ 



(19) 
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Let us define, for future convenience, the linkQ element of exp(So) as follows: 



3=0 



hj+i (A 
hiPt) 



(20) 



It is now evident that this basic element, which will be denoted also as = 
C°(9g, 9g + i), depends only on 6g, 9g + i, which are the invariant angles for the Polyakov 
lines Pg, Pg+i in the sites joined by the link. Indeed from now on we will always 
assume to have gauge-rotated the Polyakov lines to be diagonal: 



Pv 







(21) 



With these definitions the zero th -order action (|l^) is simply given by 

exp(S )=n^ • 



(22) 



However let us stress that the action that was generally used in the previous 
attempts to obtain mean field estimates of the deconfinement temperature, was 
actually a simplified version (truncated at the first representation) of eq. ([TDD : 



S p(Pt) = J2\fttJ2 cos{9g) cos(6^ H 



(23) 



It is interesting to notice that in the N t = 1 case the character expansion 
contained in eq . (|19D can be summed exactly This can be easily seen by writing the 
explicit form for the characters in eq.(|I^) : 



exp(So) = II 1 + E 



Then using the relation: 



hj+i(Pt 
hiPt) 



sin[(2j + 1)0,,] sin[(2j + l)fl gH 
sm(9g) sm(9g +i ) 



(24) 



2 sin[(2r + l)fe] sin[(2r + l)0 3+i ] = cos[(2r+l)(9 3 -9g +i )} -cos[(2r+ 1) + 

(25) 

and the well known expansion: 



J>«»° = J (/?)+ 2 £ !*(/?) cos(A;0) , 

k=l 



it is easy to obtain: 



exp(So) = II 



o Pt cos(6 s -6 s+i ) _ p p t cos{e 3 +6 3+i ) 



Ah((3 t )sm(9 £ ) sm(9 s 



(26) 



(27) 



1 The links we are referring to are those of the (i-dimensional spatial lattice, corresponding to 
one "horizontal" slice in the original d + 1-dimensional lattice 
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3.3 First order approximation 



The 0(j3g) effective action (p2|) contains just nearest-neighbour interactions between 
the Polyakov loops. We shall show below that 0(01) contributions to the effective 
action are of plaquette type, namely they involve all the four invariant angles of the 
Polyakov lines which belong to a given plaquette. This interaction is more general 
than the nearest-neighbour one, but it is still short ranged, in agreement with the 
hypothesis which is at the basis of the SY conjecture discussed above. 



3.3.1 The adjoint representation term. 

To calculate the contribution coming from the adjoint representation term, we have 
to select in eq.flTBD the term: 



X 



[28) 



To study the integral (p8|), we first note that all the spacelike plaquettes in the 
same spatial position give evidently the same contribution, regardless of the time 
t; therefore the sum over the time positions in ( p8|) simply results in a N t factor. 
Secondly, it is convenient to use the following relation for the SU(2) characters: 



Xi = (xi) 2 - 1 



(29) 



The "—1" simply reproduces the zeroth order term, and gives a renormalization 
of order fi 2 a to such contribution. The integral along the plaquette can now be 
decoupled into integrals over a single link matrix, by writing explicitly lx±(Ug )t ;ij)] 2 
as a product of elements (in the fundamental representation) of the link matrices. 
Thus eq. (^) can be rewritten in terms of the following integrals over the unitary 
spacelike link matrix U : 



B, 



a/375 



DU 



( 



\ 



hj+i(Pt 

"MAT 



Xi 



(JJP^Pl) U aP Ul s (30) 



where a, . . . = 1,2 are the indices of the U matrix in the fundamental representation. 
Making use of the invariance of the measure and of the argument of Xj m ec i- (|30| ) 
under the transformations 



U, 



a/3 



act U a p , 



and 



U a f3 ~^ U aj 3 LOpp 



u. 



7<5 



-1\ 



hi ^7<5 



(31) 
(32) 
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(where the diagonal unitary matrix uj has the property uj 2 = 1), one can conclude 
that -B Q/ 3 7 5 = unless (3 — 7 and 5 = a. As a consequence, the integral ( |3~0"D depends 
on the invariant angles of the Polyakov loop only and can be written as follows: 



B 



a(3-yS 



\Pxi Px+i) = B, 



Sa^Ssn, C. 



(33) 



(no summation over repeated indices). Moreover, it is not difficult to show that C a p 
is a real symmetric matrix. By using these results, we can write the contribution 
([2§|) to the effective action only in terms of the invariant angles of the Polyakov line 
Pg. Eq. (E3) becomes: 



3iV y 



IT Tr [C(9$, 6x+i)C(9g+i, 9g+i+j)C '{9g+i+j , 9g+j)C(9$+j, 9g)] — 1 



(34) 

The subtraction of the term (-1) in ([34] ) is due to the term (-1) in (|29|) , whereas the 
matrices C(9g, 9x + 1) are the normalized version of C: 



a 



a,f3{9g, 9g+i 



a 



a,f3\9xi 9g+i 



(35) 



The last step is the explicit evaluation of the matrix elements C a p. This calcu- 
lation is described in the Appendix. The final result is: 



C11 — C22 
C12 — C21 



2(^0 + C X ) 



(36) 



with: 



Ci 



1 



2 sin 6** sin #f 



E 



X+tj 



'x+i ^j=Q l 

+2 cos [(2j + 1)0*] cos [(2j + l)9 £+l ] 



l2j +2}JJt) 

(2j + l)/ 1 (A) iVt 

iVt-, 

+ 



'2j+l(A 



/o(A 
A (A 



(37) 



Co was defined in eq.(pO|), and, as expected, C\\ + C12 = Cq. 

Eq.([!7|) could seem a bit complicated, but it is actually very easy to implement 
it in a mean field analysis or in the SY type mapping described in the next section. 

In the N t = 1 case the sum over the representations can be performed exactly, 
and a closed expression for the C a p coefficients can be obtained. This can be done 
by using the identity: 

J(/3) n _! - I(P) n+1 = 2nI(P) n (38) 
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and eq(p6l). The result is: 

p /3t cos(e s ~e 3+t ) (3 t cos(e 3 -e 3+i ) _ p p t cos{e 3 +e 3+l ) 



Cll{ - Nt l) 4J 1 (A) sin(^) sin(fe + ,) WMPt) sm 2 (fe) sin 2 (0* +i ) 

e /3 t cos(8 s -6 s+l ) _ e l3 t cos(e s +9 s+i ) e f3 t cos(8 s +8 3+1 ) 

Cu{Nt = 1) = SfchiPt) sin 2 (^) sin 2 (^ +l ) " 4/ 1 (A) sin(^) sin(^ +i ) ' (39) 
3.3.2 Pair of fundamental representations. 

To calculate the contribution coming from a pair of fundamental representations, 
we have to select the last term in (|i~6|) : 



:0i, 



X 



x E^-Et^taXiC^tijvklC^taw) • ( 40 ) 

This term can be evaluated following the same pattern of the calculation in the 
adjoint representation case. It is however important to notice that in the N t = 1 
limit this contribution exactly disappears. Since this is the limit in which we shall 
be interested in the following we shall omit here the details of this calculation which 
can be found in 01. 



3.4 Independent approach to N t = 1 

The interesting feature of the N t = 1 limit is that in this case the theory that we 
are studying becomes a particular instance of the Migdal-Kazakov model. This 
connection was already noticed in QTUJ and was the origin of our previous analysis 
in the N — > oo limit 0. All the integrals that we have described in the previous 
sections can be directly evaluated in this particular instances of a nontrivial 

generalization of the so called Itzykson-Zuber integral, evaluated in ||11|. We refer 



the reader to |LT] for a comprehensive discussion on this interesting subject and 
simply report here the results which are useful for our analysis: 
zeroth order contribution: 



J dU 3ii exp |yT r/ (y(x)U S;i V^x + i)uV 



e Pt cos(0 3 -8 s+i ) _ e /3 t cos(8 3 +8 3+l ) 



(41) 

2Asin(^)sin(^ +J ) 1 ; 

The first order contributions can be extracted by the correlators defined in 



as: 



J dU S;i exp {f Try. (v(x)U^(x + i)U^ } 



p.o 

- ■ (42) 
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It is easy to see that these correlators must be of diagonal form, namely: 



(43) 



where the C^^{x] i) are equivalent, apart from the different normalization, to our 



Cm matrix elements. They turn out to be [11 



Ci,i(x;i) = C 2 ,2(x',i) 
Ci >2 {x;i) = C 2> i{x;i) 



2f3 t sm( 



sin 



'x+i) 



-2l3 t sin(G 3 )sin(6 3+i ) 



1 _ e -2Asin(0 2 ) s in(f? a+i )j ( 2 /5 t sin(fe) sin(0a +i )) 
x _ e -2/3 tS m(e 2 ) S in(% +i ) ^ + 2 /3 t s in(fe) sin(% +i )) 
(l _ e -2/3 tSi n(e :e )sin(e 2+l )^ (2^ s in(^) sin(fe +l )) 



.(44) 



It is now only matter of straightforward algebra (one must also take into account 
eq.(|35D and the normalization constant G{(3) defined in ( fl~4|) ) to show that these 
expressions eq. (^Tf) and ( [44]) are exactly equivalent to our results (|27|) (p9|) . 



4 Determination of T c 

As we discussed in the introduction, the standard approach to the determination 
of the critical temperature would be at this point a mean field analysis of the 
above constructed effective action. However this approach is rather unsatisfactory. 
For instance if we take the standard mean field approximation of the zeroth order 
action, truncated at the first representation, eq.fl2~3|) ( for which the calculation 
can be performed exactly, see for instance 0) it is easy to see that the resulting 
critical coupling in the N t = 1 case is j3 c = 2/d. In (3+1) dimensions we know that 
(3 C = 0.8730(2) M, and it is clear that the standard mean field derivation which 
gives in this case (3 C = 0.666... is largely unsatisfactory. It is possible to improve 
this result keeping the full effective action instead of its truncated version, and 
improving the mean field approximation. This gives a much better result, which 
however always remains between 5 and 10 % below the Montecarlo result (see M 
for a discussion of this approach). 

In this section we want to discuss a completely different approach, which makes 
explicit use of the mapping between the SU(2) model and the d- dimensional Ising 
model and turns out to be much more powerful of the mean field approach. Up 
to our knowledge the only attempt along this line was made by J. Polonyi and 
K.Szlachanyi in 0, but since they were constrained to keep in the various stages of 
their analysis only the very first order in f3 t (and to neglect (3 S ) their result was rather 
unsatisfactory. We shall review their approach below. The novelty of our approach 
with respect to this previous attempt is twofold. First, we keep all the orders in the 
(3 t expansion of the interaction. Second, and more important, we use the explicit 
knowledge of the first order in the f3 s expansion to map the original gauge model 
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to the equivalent Ising model, explicitly relating the gauge coupling and the Ising 
coupling in a new original way, completely different from that of ||. From the 
knowledge of the location of the Ising phase transition we can thus reconstruct the 
exact critical deconfinement temperature. 

Let us follow this procedure in the N t = 1 case. The crucial point of the whole 
approach is the identification of the Ising variable "embedded" in the Polyakov 
loops. Let us follow as a first example the analysis of ref. ||. The simplest way to 
extract these Ising variables is to decompose the Polyakov loops: 

P{x) = 2COS(0*) , 6jt £ [— 7T, 71*] . (45) 

which are SU(2) variables, into the product of a Z2 variable <r(x), which is simply 
a sign, and a SO (3) variable: 

P(x) = 2cos(0*) , B s e [-tt/2, tt/2] . (46) 

Then by integrating over P one ends up with the desired effective action of the 
Ising type. All these steps can be easily performed if we keep at each stage only 
the very first order in (3 t (and neglect j3 s ). In this case the first step gives us the 
effective action S p of eq.fl2"3"|). Then, it is possible to see that at the first order in 
(3 t the SO (3) variables decouple and can be integrated out exactly. The effective 
action becomes (neglecting an irrelevant overall constant): 

SM = E{/5 /sm9 E <r(S+i)\ , (47) 

X I 1=1 J 

with /3 Isin9 = ^sflt- This is exactly the ordinary Ising action which is known to 
have, both in two and three dimensions a second order phase transition located at 
Rising = Q.44068679 ■ • ■ and p 1 ^ 119 = 0.221652(3) |TJ in two and three dimensions 
respectively. From this we obtain the following values for the deconfinement tem- 
perature: fa = ^fc 3 ™ 9 , namely (3 t = 2.447 and (3 t = 1.231 in (2+1) and (3+1) 
dimensions respectively. Even if the order of magnitude is essentially correct it is 
easy to see that these estimates are even worse than the plain mean field analysis 
on the action eq.(E3|). Notwithstanding this, let us stress again that this approach 



is very interesting in itself, because it allows to have a deeper physical insight on 
the mechanism underlying the deconfinement transition. Let us now improve this 
analysis by keeping in the various step all the orders in the (3t expansion and the 
first order in f3 s . 

First, let us start from eq. (f2~7|) , which is the exact effective action, taking into 
account all orders in f3 t and let us apply the same recipe as above to extract the 
embedded Ising variables. We obtain after some simple algebra: 

(d ( e Pt sin(%) sin(0 a+i ) _ g -/3t s'm(e s ) sin(6» 2+ j) ~j \ 
n( — w^mmim — s) 

x e f3 t cos(e s )cos(8 s+l )<T(x)cr(x+i) ^gN 
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Let us define 

a(x, i) = cos(6g) cos(^ + i) , 

j(x, i) = sm(6g) sin(6» £+i ) . 

Since at the deconfinement point the SO (3) variables are not critical, we can 
assume that 8(x) is described also at the deconfinement point by a constant master 
field 9q and that fluctuations in 6 can be neglected. Let us consistently define: 



The problem is thus reduced to the identification of this master field, in terms 
of which we could write f3 t)C = (3 I c sm9 / olq 

The only way we have to find this master field is to identify also at the level 
of the O(Pg) effective action the underlying Ising model. However this is clearly a 
non-trivial task. In fact at a first glance one could think that it is impossible to 
map our 0(/3f) effective action into an equivalent Ising model because the O(0f) 
contribution is not of the nearest-neighbour type. For instance, if we separate also 
in this case the Z 2 degrees of freedom from the SO (3) ones as we did above, due 
to the plaquette structure of the action, the Z 2 degrees of freedom exactly cancel 
each other, and the underlying Ising model seems to be lost. However there is a 
completely different way to recognize such an Ising model in the 0(/3 2 ) interaction. 
In fact, a remarkable and non-trivial consequence of the particular form of the 0((3f) 
interaction is that it can be exactly reorganized as the first term of the strong 
coupling expansion of an ordinary nearest-neighbour Ising model, thus allowing 
to complete the identification. In fact, if we keep as above a constant value of 
9 = 9 , then the relations = C 2j 2 and Ci j2 = C 2i i, and above all the fact that 
Ci 5 i + Ci j2 = 1, allow us to interpret eq. (|3~4D as the strong coupling expansion of 
a (^-dimensional Ising model with action, say, Sj S i ng = J <j(x)<j(x + i) where the 
values i = 1 and i = 2 of the indices of C^j denote the +1 and —1 values of the 
Ising spin a(x) and the Ising coupling J is given by 



This induced Ising model can be considered as the replica at the first order in 
(3 S of that described above at the zeroth order in f3 s . Again we must require the 
coupling J, to be at its critical value J c . Solving eq.(|49"D with respect to 70, we find 
70 = 1.3957/ /3 £)C in (2+1) dimensions and 70 = 0.67383//5 tjC in (3+1) dimensions. 
Combining these values with the above discussed relation: (3 t)C = Pc Szn9 /o!o, we 
finally find: 

(2+1) dimensions: 



a = cos 



.2 



70 = sin 2 (6>o 




(49) 



(3 t = 1.836 , a = 0.2398 



0.337tt 
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(3+1) dimensions: 

Pt = 0.8954 , a = 0.24 75 , 9 = 0.334tt . 

As we anticipated above, these values for f} c are much lower than those of ref. 0, 
and the one for d = 3 is in good agreement with the result obtained by Montecarlo 
simulations. 

5 Conclusions 

The approach outlined in sect. 4 can be extended also to N t > 1 |J. In following 
this extension one must take care of some non-trivial features of the models, like 
the fact that the critical coupling (3 C as a function of N t obeys different scaling laws 
in (2+1) and (3+1) dimensions. However the pattern of our approach needs not to 
be changed. The agreement with the Montecarlo results (when they exist) remains 
very good. Since our approach is not limited by the magnitude of N t , we can hope 
that, as N t increases, a sensible continuum limit for T c could be taken. To reach this 
goal we must be able to reconstruct the correct scaling laws in the large N t limit. 
This is certainly possible for the (2+1) dimensional model (see || for details), but 
it is still an open problem in the (3+1) dimensional case. In any case, besides the 
numerical results, we think that the improvements that we have discussed in this 
contribution both in constructing the effective action and in extracting the critical 
coupling can help us to have a better and deeper understanding of the physics of 
the deconfmement transition in lattice gauge theories. 
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Appendix 

In this Appendix we shall evaluate the matrix elements Cm- 
To this end let us select in the sum over representations contained in Cm the j th 
term: 

A j kl = J D U\U kl \ 2 Xj (V(x) U V\x + i)U^) (A.l) 
so that Cm can be written as: 

N t 

A r k i ■ (A.2) 



a 



kl 
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Let us use the following relation: 



(where [j] denotes the integer part of j), which is a direct consequence of the 
identification of the SU(2) characters with the Tchebichef polynomials of second 
kind: Xj(U) = U 2 j(cos(9)) (where 9 denotes, as usual, the invariant angle of the 
matrix U). 

We can rewrite ( |A.1| ) as: 



ki - S / D U\U^^(V{x) U V*@ + i)tf) . (A.4) 

Since the U matrix elements always appear in the form \U kl \ 2 (with the indices 
in the fundamental representation) it turns out that a very useful parametrization 
is: 

3 

U = a l + ^aj(Tj (A. 5) 

where Oi are the Pauli matrices, the are real numbers constrained by: J2i=o a 2 = 1. 
In this parametrization we have |f/ 12 | 2 = a\ + a\ and \U n \ 2 = \U 22 \ 2 = a 2 , + a\. 
Setting a\ + a 2 . = x we see that the xi in eq. (|A.4j) becomes 

Xi(V(x)UV*(x + i)rf)=g + hx (A.6) 

with g = 2cos(6^ — fe+j) and h = —Asm8g sin^ + j. The measure DU becomes dx, 
with integration limits and 1, according to the above mentioned constraint on the 
a«. The A 3 kl integrals become: 

U] (-l) fc (2j-A;)! 



Before evaluating these integrals, as a preliminary exercise, let us calculate the 
simpler integral in which no contribution coming from the spacelike plaquette is 
present. Let us call it A 3 Q : 



If we are able to evaluate the integrals and sum up the series we shall find an 
alternative way to go from eq. (|i"8f) to eq.(|K}). 



16 



The integrals in the sum of eq.( [A.9| ) can be done straightforwardly: 

(g + h) n+l - g n+1 



dx(g + hx) 1 



h(n + 1) 



(A.10) 



inserting this result in eq. (|A.9| ), and using the explicit expression for g and h we 
have: 



■ = |l (-l) fc (2j - fc)! [2 cos(fe + 9 s+l )] n+l - [2 cosCflg - 9 s+i )] n+1 
t k\(2j + l-2k)\ 



—4 sin Qg sin 6. 

By using the explicit expression of the Tchebichef polynomials of first type: 

V2] f_ 1 )k (n _ k _ 1)l 



(A.11) 



T n (cos(9)) = cos(n9) = 77 H 



k=0 



k\ (n-2k)\ 



[cos(n6)} 



2j-2k 



(A.12) 



we can rewrite A J as 



A 1 
^0 



cos((2j + l){9 3 - 3+i )) - cos((2j + l)(0g + e 3+i )) 



2(2j + l)sin^ sin^ 



(A.13) 



x+i 



using eq.(p5l), inserting the result in the sum on the representations, and using the 
explicit expression for the characters Xj we exactly obtain eq. flT9|) . 

Let us now calculate A\ x . Also in this case the integral contained in eq. (|A.7|) 
are straightforward: 



dx(l - x)(g + hx) n 



(g + h) n+2 - g n+2 



(A. 14) 



h 2 (n + l)(n + 2) h{n + l) 

The second term in eq.( |A.14|) can be treated exactly as we did above for A J Q . 
The contribution to A 3 U coming from it is (see eq. ( |A.13|) ): 

cos((2j + l)(^-^ +J )) 



2(2 j + 1) sm6g sin 9g +i 
The first term of eq.( |A.14j) , after inserting the expression for g and h gives: 



(A.15) 



b1 

£ kl (2J + 2 



l) h (2j - k)\ [2 cos(fe + fe +4 )]"+ 2 - [2 cos(fe - 



n+2 



2A;)! 



16(sin6 l £ ) 2 (sin6'^ 



(A.16) 



To sum this series the following trick is needed. Let us divide and multiply for 
2j + 1 and let us split this factor at the numerator as (2j + 1 — k) + k. Then the 
sum ( [A.16 ) is splitted into two sums that, after suitable redefinition of the indices 
can be reduced to the sum ( |A.12| ). Collecting together all the pieces one finds: 
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M 



cos((2j + 2){6 S + 9 iM )) - cos((2j + 2){6 S - 9g +i )) 



11 



8(2j + l)(2j + 2)(sinfe) 2 (sin# 



X+lJ 



+ 



cos((2j)(^ + 9 g+i )) - cos((2j)(^ - Og+i)) 

8(2j + l)(2j)(sin^) 2 (sin^ +i ) 2 
co8({2j + 1)(9 S - 6 g+i )) 
2(2j + l)sin0 a sinfl 



(A. 17) 



x+i 



By using the definitions ( A.7 ),( A~lj|) and ( |A.9|) one can immediately obtain A\ 2 
as the difference: A\ 2 = A J — A J n . 

Inserting these results in eq. (|A.2|) one finally obtains the C^i coefficients. Simple 
trigonometric relations allow to write these coefficients in a compact form: 



Cii — C 2 2 — + 



(A.18) 
(A.19) 



with: 

Ci 



1 



2 sin 9? sin 9 



Xj(9z)Xj(9g+i 



<-2j+2\fit 



hj(Pt 



(2j + l)/ 1 (/? t )^ 



+2 cos [{2j + 1)0*] cos [(2j + 1)0**] 



hj+i (fit 
, hiPt) 



+ 



1 N t 



A.20) 
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